Spontaneous and induced dynamic fluctuations in glass-formers I: 
General results and dependence on ensemble and dynamics 
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We study theoretically and numerically a family of multi-point dynamic susceptibilities that quan- 
tify the strength and characteristic lengthscales of dynamic heterogeneities in glass-forming mate- 
rials. We use general theoretical arguments (fluctuation-dissipation relations and symmetries of 
relevant dynamical field theories) to relate the sensitivity of averaged two-time correlators to tem- 
perature and density to spontaneous fluctuations of the local dynamics. Our theoretical results are 
then compared to molecular dynamics simulations of the Newtonian, Brownian and Monte-Carlo 
dynamics of two representative glass-forming liquids, a fragile binary Lennard- Jones mixture and a 
model for the strong glass-former silica. We justify in detail the claim made in [Science, 310, f 797 
(2005)], that the temperature dependence of correlation functions allows one to extract useful infor- 
mation on dynamic lengthscales in glassy systems. We also discuss some subtle issues associated to 
the choice of microscopic dynamics and of statistical ensemble through conserved quantities, which 
are found to play an important role in determining dynamic correlations. 

PACS numbers: 64.70.Pf, 05.20.Jj 



I H S [2i|, and experi- 
2j,|25|,l26|,|27|, of various aspects of dynamic 
heterogeneity has greatly advanced our understanding of 
the behavior of systems close to the glass transition. In 
particular, multi-point susceptibilities have been devised 
to quantify the behavior and magnitude of the putative 
growin g dy namical lengthscale [1, [H, [H, [H, [l^, [H, 
H, [30I. IsTl [3^ . m, [U, and experimental studies have, 
for several materials, directly determined the number of 
molecular units that move coo per atively near the glass 
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Diverse materials, ranging from molten mixtures of 
metallic atoms, molecular and polymeric liquids, and col- 
loidal suspensions may form glasses if sufficient under- 
cooling or densification occurs P, 0, Q . A glass may be 
characterized mechanically as a solid, but microscopically 
lacks the long-range order of a crystal. Close to vitrifi- 
cation, the viscosity of glass-forming systems increases 
dramatically and sensitively as the thermodynamic con- 
trol variables are changed. Furthermore, some degree 
of universality is observed in the thermal and temporal 
behavior of systems close to the glass transition, even 
though the material properties of such systems may be 
vastly different [l|, 0]. Despite decades of intense theo- 
retical and experimental work, the underlying causes of 
this interesting behavior are not well understood. 

The observed quasi-universal behavior of glassy sys- 
tems might be related to the existence of a growing 
lengthscale as the glass transition is approached. The 
search for such a correlation lengthscale has led to in- 
tense activity in recent years. Static structural indicators 
have repeatedly failed to show any evidence of collective 
behavior. Indeed, the static structure of a supercooled 
liquid hardly differs from that of the same liquid at rel- 
atively high temperatures. Clearly all simple structural 
correlations remain short-ranged as the glass transition 
is approached 3- It has become manifest in the last 
decade that interesting behavior is revealed by spatially 
correlated dynamics. As a whokjSuch effects are referred 
to as 'dynamical heterogeneity' [1, H, 0, H, 111 ■ 

The investigation, via theory [13, [O, [H, [H, [H, 



Despite recent 



transition [22, |23|, |25|, |26|, |3£ 

breakthroughs, much more work needs to be done to 
fully characterize such a behavior both experimentally 
and theoretically. 

In the present work, contained here and in a compan- 
ion paper [3^, we make a step towards this goal by in- 
vestigating in detail different susceptibilities that may 
be categorized according to the induced or spontaneous 
nature of the measured fluctuations [37{ . Spontaneous 
dynamic fluctuations can be characterized by four-point 



functions, as prop osed and studied earlier [a [H, [If 
[20,M|2a[33,[3l[[33,[3l,[3i|, or three-point functions, as 
in [371. |39| . Instead, fluctuations can be induced hy mon- 
itoring the change of dynamical correlators that follows 
a change of an external control parameter, e.g. temper- 
ature [33, H^. As we shall show, it is possible to re- 
late induced and spontaneous dynamical fluctuations via 
fluctuation-dissipation relations as proposed in |33| • This 
provides a very valuable experimental tool to measure dy- 
namic fluctuations since, as usual, induced fluctuations 
are much easier to measure than spontaneous ones. 

Using molecular dynamics simulations of different 
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archetypal glass-forming liquids (e.g. "strong" materials 
that exhibit an Arrhenius temperature dependence of the 
viscosity, and "fragile" ones, whose viscosity displays a 
super- Arrhenius temperature dependence) we shall show 
that in the slow dynamical regime a considerable fraction 
of spontaneous fluctuations can be attributed to energy 
fluctuations: since the dynamics is spatially correlated a 
local energy fluctuations induces a change in the dynam- 
ics over a much larger range. 

Our analysis will however reveal that global four-point 
correlations describing the fluctuations of intensive dy- 
namical correlators may depend both on the statisti- 
cal ensemble and on the underlying microscopic dynam- 
ics. Local four-point correlations measuring the corre- 
lation between the relaxation dynamics at a finite dis- 
tance apart of course do not depend on the statistical 
ensemble but they do depend on the underlying micro- 
scopic dynamics. This is striking because it is known 
that correlators measuring the average dynamics do not 
depend in the relevant glassy regime on the microscopic 
dynamics (Newtonian or stochastic [iO, |4l| ) . We address 
this problem both theoretically and numerically, and con- 
clude that, although the underlying physical mechanisms 
are the same, dynamical correlations depend quantita- 
tively on the conserved physical quantities (and global 
four-point correlations even on the statistical ensemble). 
For example, the absolute magnitude of global sponta- 
neous dynamical fluctuations in a Lennard-Jones system 
in the NVT ensemble obtained from Brownian dynam- 
ics (BD) or Monte-Carlo dynamics (MC) are very simi- 
lar, but are considerably smaller than that obtained with 
Newtonian dynamics (ND) in the same NVT ensemble, 
whereas ND simulations performed in the NVE ensem- 
ble yield results that are close to the BD and MC results 
in the NVT ensemble. However, we stress that all our 
results point toward the conclusion that the behaviour of 
all these quantities as the glass transition is approached is 
governed by the growth of a unique dynamic correlation 
length, at least in the numerically accessible regime. 

The aim of the present paper is to provide the reader 
with the physical picture underlying the dynamical sus- 
ceptibilities introduced in (37| , along with more technical 
elements based on general field-theoretical considerations 
and detailed numerical investigations of different real- 
istic glass- forming liquids. In a companion paper [ssj . 
we present some quantitative predictions for these sus- 
ceptibilities, obtained within different theoretical mod- 
els: mean- field spin glass models [i^l, mode-coupling the- 
ory [4^ , and kinetically constrained models [i^l which we 
again confront with the results from molecular dynamics 
simulations. The present paper is arranged in four sec- 
tions. In Section |TT] we present the physical motivations, 
definitions and physical content of several multi-point dy- 
namic susceptibilities. We derive in particular general 
results for the ensemble dependence of dynamic fluctu- 
ations, fluctuation-dissipation relations, and bounds be- 
tween induced and spontaneous dynamic fluctuations. In 
Section [ni] we present a field-theoretic derivation of the 



behavior of dynamic fluctuations for various types of mi- 
croscopic dynamics. This is particularly useful in identi- 
fying the precise physical mechanism leading to a growth 
of dynamic correlations, and the dependence of multi- 
point susceptibilities on the microscopic dynamics. In 
Section [IVI we summarize our various theoretical predic- 
tions and extract some important consequences, relevant 
to experiments, that need to be tested numerically. In 
Section |V] we present the results of detailed molecular 
dynamics simulations of two model glass-forming liquids, 
a fragile binary Lennard-Jones mixture and the strong 
BKS model for silica. We compare spontaneous and in- 
duced fluctuations and show that, as predicted theoreti- 
cally, dynamic correlations strongly depend on the choice 
of microscopic dynamics and statistical ensemble. Our 
results suggest, however, that a unique dynamical length- 
scale governs the growth of dynamical susceptibilities in 
all cases. In Section IVII we give the conclusions of our 
study. Although very natural in spin-systems, four-point 
correlators in liquids mix dynamical heterogeneities with 
different physical effects (in particular, energy and den- 
sity conservation) and might therefore not be the most 
effective object to work with. On the other hand, we fully 
confirm the claim made in (37j , that the temperature de- 
pendence of correlation functions allows one to extract 
rich and useful information on dynamic lengthscales in 
glassy systems [45 1. 



II. 



MULTI-POINT DYNAMIC CORRELATORS 
AND NEW LINEAR SUSCEPTIBILITIES 



A. Why four-point correlators? The spin glass case 

No static correlation has yet been found to reveal an 
notable feature upon approaching the glass transition 
licj . Any lengthscale associated with the slowing down 
of the system must therefore be hidden in some dynamic 
correlation function. This issue is in fact deeply related 
to one of the most important question pertaining to the 
physics of disordered systems - how can one define long- 
range amorphous order in such systems? 

We know from the theory of spin glasses that the above 
oxymoron has in fact a precise answer: some hidden 
long-range order indeed develops at the spin glass tran- 
sition [47[. In order to reveal this long-range order, con- 
ventional two-point functions are useless. Even if spins 
Sx and Sx+y have non-zero static correlations (sx^x+y) 
in the spin glass phase, the average over space for a given 
distance y vanishes because the pairwisc correlations ran- 
domly change sign whenever x changes. The insight of 
Edwards and Anderson is that one should first square 
(sxSx+y) before averaging over space [l^. In this case, 
the resulting (four-spin) correlation function indeed de- 
velops long-range tails in the spin glass phase. This cor- 
relation in fact decays so slowly that its volume integral, 
related to the non-linear magnetic susceptibility of the 
material, diverges in the whole spin glass phase [49j . 
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The Edwards- Anderson idea can in fact be understood 
from a dynamical point of view, which is important for 
understanding both the physics of the spin glass just 
above the transition, and its generalization to structural 
glasses. Consider, in the language of spins, the following 
four-point correlation function: 

Siiy, t) = [{S^{t = 0)s^+y{t = 0)SxWSx+y(i)>]x, (1) 

where the brackets [...]x indicate a spatial average. Sup- 
pose that spins Sx and s^+y develop static correlations 
(sxSx+y) within the glass phase. In this case, 6*4(7, i ^ 
00) will clearly converge to the spin glass correlation 
[{sxSx+yy^]x- More generally, 5*4 (y,t) for finite t is able 
to detect transient tendencies to spin glass order, for ex- 
ample slightly above the spin glass transition tempera- 
ture Tc- Close to the spin glass transition, both the per- 
sistence time and the dynamic length diverge in a critical 
way: 

S4y,t)^y^-^-^s(^^,'-y (2) 

where ^ ^ (T - Tc)-" and r - (T - T^)-'-'' . As men- 
tioned above, the static non-linear susceptibility diverges 
as / dySi{y,t 00) ~ More generally, one can 

define a time-dependent dynamic susceptibility as: 

Xi{t) = j dy S4{y,t), (3) 

which defines, provided 6*4(0, t) ~ 1, a correlation vol- 
ume, i.e. the typical number of spins correlated in dy- 
namic events taking place over the time scale t. As we 
shall discuss below, Xii'^) can also be interpreted as a 
quantitative measure of the dynamic fluctuations. Note 
however that the precise relation between Xi and ^ de- 
pends on the value of the exponent 77, which is physically 
controlled by the detailed spatial structure of S4: 

X4(i = r) « e'""- (4) 

Therefore, spin glasses offer a precise example of a sys- 
tem which gets slower and slower upon approaching Tc 
but without any detectable long-range order appearing 
in two-point correlation functions. Only more compli- 
cated four-point functions are sensitive to the genuine 
amorphous long-range order that sets in at Tc and give 
non-trivial information even above Tc- In the case of spin 
glasses it is well established that the transition is related 
to the emergence of a low temperature spin glass phase. 
In the case of the glass transition of viscous liquids the sit- 
uation is much less clear. There might be no true phase 
transition toward a low temperature amorphous phase. 
It is still reasonable to expect that the dramatic increase 
of the relaxation time is due to a transient amorphous or- 
der that sets in and whose range grows approaching the 
glass transition. Growing timescales should be somehow 
related to growing lengthscales [13 . A good candidate 



to unveil the existence of this phenomenon is the func- 
tion 5*4(7, introduced previously, since nothing in the 
above arguments was specific to systems with quenched 
disorder. The only difference is that although transient 
order is detected in Si{y,t) or its volume integral Xi{t) 
for times of the order of the relaxation time, in the long 
time limit these two functions may not, and indeed do 
not in the case of supercooled liquids, show long-range 
amorphous order. This roots back to the different nature 
of the glass and spin glass transitions (see the discussion 
in 



B. Supercooled liquids and more multi-point 
correlations 

In the case of liquids, we may consider a certain space 
dependent observable o(x, f), such as, for example, the 
local excess density (5p(x, i) ~ p{'X-,t) — p^, where po is 
the average density of the liquid, or the local dipole mo- 
ment, the excess energy, etc. We will assume in the fol- 
lowing that the average of o(x, t) is equal to zero, and 
the variance of o(x, t) normalized to unity. The dynamic 
two-point correlation is defined as: 

Co(r,t) = [o(x,t = 0)o(x + r,t)]„ (5) 

where the normalization ensures that Co(r = 0, t = 0) = 
1. The Fourier transform of Co(r, t) defines a generalized 
dynamic structure factor S*o(k, t) [5l|. All experimental 
and numerical results known to date suggest that as the 
glass transition is approached, no spatial anomaly of any 
kind appears in Co(r, t) (or in S*o(k, t)) although of course 
there could still be some signal which is perhaps too small 
to be measurable. The only remarkable feature is that 
the slowing down of the two-point correlation functions 
often obeys, to a good approximation, "time-temperature 
superposition" in the a-relaxation regime t ^ Ta, i.e.: 

Co(r,t)«q„(r)/(^-A_^, (6) 

where qo is often called the non-ergodicity (or Edwards- 
Anderson) parameter, and the scaling function /(x) de- 
pends only weakly on temperature. This property will be 
used to simplify the following discussions, but it is not a 
crucial ingredient. 

Whereas Co(r, t) measures how, on average, the dy- 
namics decorrclates the observable o(x, i), it is natural 
to ask whether this decorrelation process is homogeneous 
is space and in time. Can the correlation last much 
longer than average? In other words, what is the dis- 
tribution (over possible dynamical histories) of the cor- 
relation Co(r,i)? Clearly, since Co(r,t) is defined as an 
average over some large volume V, the variance of 
C*o(r, t) is expected to be of order /V, where ^ is the 
Icngthscalc over which Co(r, t) is significantly correlated. 
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More precisely we define: 

= I Y^o{^,0)o{^ + r,t)oi^',0)oi^' + r,t) 



(7) 



which, using translational invariance. can be transformed 
into the space integral of a four-point correlation: 



(8) 



where 



54(y, t) = <j [o(x, 0)o(x + r, t)o(x + y, 0)o(x + y + r, t)]^ 
-[o(x,i = 0)o(x + r,<)]2l (9) 



The variance of Co(r, t) can thus be expressed as an inte- 
gral over space of a four-point correlation function, which 
measures the spatial correlation of the temporal correla- 
tion. This integral over space is also the Fourier trans- 
form of 5*4 (y,i) with respect to y at the wavevector q 
equal to zero. We want to insist at this stage that r 
and y in the above equations play very different roles: 
the former enters the very definition of the correlator we 
are interested in Eq. ([5]), whereas the latter is associated 
with the scale over which the dynamics is potentially cor- 
related. Correspondingly, great care will be devoted in 
the following to distinguish the wavevector k, conjugate 
to r, and q conjugate to y. 

Specializing to the case r = (local dynamics), one 
finally obtains [13]: 



y2 



N ' 



(10) 



The analogy with spin glasses developed above sug- 
gests that this quantity reveals the emergence of amor- 
phous long-range order; it is in fact the natural di- 
verging susceptibility in the context of p-spin descrip- 
tions of supercooled liquids, where a true dynamical 
phase transition occurs at a certain critical tempera- 
ture [13, m, m, m, \E^. since in real systems no true 
phase transition is observed, one expects Xii^) to grow 
until t ~ Ta and decay back to zero thereafter. Until 
Tc, there cannot be strong differences between a system 
with quenched disorder and a system where disorder is 
dynamically self-induced. 

However, contrary to spin glasses, for which an un- 
derlying lattice structure exists, viscous liquids consist 
of molecules or atoms having continuum positions. As 
a consequence, one has to coarse-grain space in order to 
measure the fluctuations of the local relaxation dynamics. 
Local now means on a region of the order of the interpar- 
ticle distance. Therefore, generically, Xii^) ~ corre- 
spond either to the fluctuations of the Fourier transform 
of Co(r, t) evaluated at a wave- vector, ko, of the order of 



the first peak in the structure factor [3^, or to a spatial 
average / drCo{r,t)w{r) where w{r) is an overlap func- 
tion equal to one for lengths of the order of 27r/fco and 
zero otherwise Q. The dependence of dynamical corre- 
lations on the coarse-graining length has been recently 
studied in [s^l and is also discussed in the companion 
paper fssj . 

Although readily accessible in numerical simulations, 
is in general very small and impossible to measure 
directly in experiments, except when the range of the 
dynamic correlation is macroscopic, as in granular mate- 
rials [36j or in soft glassy materials where it can reach the 
micrometer and even millimeter range [ssl . [53 | . The cen- 
tral idea of this work is that induced dynamic fluctuations 
are more easily accessible than spontaneous ones, and 
can be related to one another by fluctuation-dissipation 
theorems. The physical motivation is that while four- 
point correlations offer a direct probe of the dynamic 
heterogeneities, other multi-point correlation functions 
give very useful information about the microscopic mech- 
anisms leading to these heterogeneities. For example, one 
expects that the slow part of a local enthalpy (or energy, 
density) fluctuation per unit volume (5ft, at x and time 
t = triggers or eases the dynamics in its surroundings, 
leading to a systematic correlation between (5/i(x, t ^ 0) 
and o(x',t = 0)o(x' + r, Tq,). This defines a family of 
three-point correlation functions that relate thermody- 
namic or structural fluctuations to dynamics. Interest- 
ingly, some of these three-point correlations are both ex- 
perimentally accessible and give bounds or approxima- 
tions to the four-point dynamic correlations. The reason 
is as follows. In the same way that the space integral 
of the four-point correlation function is the variance of 
the two-point correlation, the space integral of the above 
three-point correlation is the covariance of the dynamic 
correlation with the energy fluctuations [sll : 



-'CH 



J dx<ix'o(x' + r,t)o(x',0)(5/i(x,0) 

^ jfj rfyKx + y + r,0o(x + y,0)<5/i(x,0)],. (11) 

Hence, using the fact that the enthalpy fluctuations per 
particle are of order y/cpkBT (where cp is the specific 
heat in units), the quantity NY^c'H / ^/cpksT defines 
the number of particles over which enthalpy and dynam- 
ics are correlated. Of course, analogous identities can 
be derived for the covariance with density (and energy) 
fluctuations. 

Now, on very general grounds, the covariance obeys 
the Cauchy-Schwarz bound: 'S^j^ < EpS^^, where 
is the variance of the enthalpy fluctuations, equal to 
cp{kBTf/N in the NPT ensemble, N = poV being the 
total number of particles. Therefore, the dynamic sus- 
ceptibility X4(0 is bounded from below by: 



X4(t) = A^E 



> 



A^E, 



CH 



cpiksT) ^ 



(12) 



5 



where, as we show below, the right hand side can be ac- 
cessed experimentally. We then discuss in Sec. Ill El how 
the above bound can be interpreted as an approximation, 
with corrections that can be physically estimated. Note 
that we chose here Xii^) to define a number of particles; 
one can of course convert it into a volume by multiplying 
Xi by vq = 1/poj the average volume per particle. Note 
also that here and in the following we will work in the 
NPT ensemble which is the relevant ensemble for exper- 
iments on molecular liquids. We will discuss later the 
generalization to different ensembles. 



C. A dynamic fluctuation-dissipation theorem and 
growing lengthscales 

Consider a system in the grand-canonical NPT en- 
semble. The probability of a given configuration C is 
given by the Boltzmann weight exp(— /3i7[C])/Z, where 
P = l/ksT and Z is the grand-partition function. Sup- 
pose one studies an observable O with the following prop- 
erties: (i) O only depends on the current microscopic 
configuration C of the system and (ii) O can be written 
as a sum of local contributions: 



O 



dx ofx) 



(13) 



In this case, a well-known static fluctuation-dissipation 
theorem holds ISlI: 



dx (o(x)5/i(0)) 



'NY.. 



OH, 



(14) 



where we decomposed the enthalpy in a sum of local con- 
tributions as well [5l| . 

Interestingly, in the case of deterministic Hamiltonian 
dynamics, the value of any local observable o(x, t) is in 
fact a highly complicated function of the initial configu- 
ration at time t ~ 0. Therefore, the correlation function, 
now averaged over both space and initial conditions can 
be written as a thcrmodynamical average: 

Co(r, t; T) = -^ly^ J dx o(x + r, i)o(x, t = 0) 



X exp 



-P J dx' h{x',t = 0) 



(15) 



Hence, the derivative of the correlation with respect 
to temperature (at fixed volume) directly leads, in the 
case of purely conservative Hamiltonian dynamics, to 
the covariance between initial energy fluctuations and 
the dynamical correlation. Defining 5't(x, t) — (o(x -I- 
r,t)o(x,0)J/i(0,0)), one finds: 



dCo{r,t;T) 1 



dT 



dxST{^,t)=XTir,t). (16) 



Hence, the sensitivity of the dynamics to temperature 
XT is directly related to a dynamic correlation. This last 



equality, although in a sense trivial, is one of the central 
result of this work. It has an immediate deep physical 
consequence, which is the growth of a dynamical length 
upon cooling in glassy systems, as we show now. 

Define Ta{T) such that Co(0,t = Ta;T) = e^^ (say). 
Differentiating this definition with respect to T gives 

^_dT^ dCoiO,t = T„-T) ^ dCo{0,t ^ To,;T) ^^^^ 



dT 



dt 



dT 



Since Co(0,t;T) decays from 1 to zero over a time scale 
Tq., one finds that generically, using Eq. (|16p : 



rfx 



(o(x,t = Ta)o(x,0)(5/i(0,0)) 
po^/c^kBT 



1 d In T„ 



Pay/cp dlnT 

(18) 

Now, Sh is of order pQy^cpksT and (o^) is normalized to 
unity, and the quantity xo = STiO,Ta)/poy^cpkBT can- 
not appreciably exceed unity. The above integral can be 
written as Xo'^^t-, which defines a volume vt over which 
enthalpy fluctuations and dynamics are appreciably cor- 
related. Note that the interpretation of vt as a true 
correlation volume requires that xo be of order one, and 
its increase is only significant if xo is essentially tempera- 
ture independent. If this is not the case, then the integral 
defined in (18) could grow due to a growing xo and not a 
growing length, which would obviate the notion that vt 
is a correlation volume. For now wc will assume these 
properties hold, and return to this crucial point theoret- 
ically in more detail in Sec. IV and with direct numerical 
evidence in Sec.V. 

Assuming Xo !i 1; a divergence of the right hand side 
of the equality (|16p necessarily requires the growth of vt- 
More precisely, as soon as increases faster than any in- 
verse power of temperature, the slowing down of a Hamil- 
tonian system is necessarily accompanied by the growth 
of a dynamic correlation length. However, as already 
mentioned above, the precise relation between dC/d\nT 
and an actual lengthscale, ^, depends on the value and 
structure of the spatial correlation function (for example 
the value of xo and the exponent 77). In the simplest case 
of an exponentially decaying S't(x, t), one finds: 



^ dCo{r,t;T) 3 
T = Stt^xoPqC 



(19) 



It is instructive to study the case of a strong glass- 
former, for which the slowing down is purely Arrhenius, 
i.e Tq = Tq exp[A /{ksT)], where A is some activation 
barrier. The volume vt is then given by: 



Vt 



d In Tn 



dhiT 



ksT' 



(20) 



which increases as the temperature is decreased, and di- 
verges as T ^ 0. This is at first sight contrary to intu- 
ition since simple barrier activation seems to be a purely 
local process. However, one should remember that the 
dynamics strictly conserves energy, so that the energy 
used to cross a barrier must be released from other parts 
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of the system. This release necessarily induces dynamic 
correlations between the Arrheniusly relaxing objects. 
We conclude that even in a strong glass the Arrhenius 
slowing down is necessarily accompanied by the growth 
of a dynamic lengthscale. Note again that this conclu- 
sion relies on subsidiary conditions that must be met. 
Indeed, it is not difficult to find examples of model New- 
tonian systems for which \dlaTa/dhiT\ grows substan- 
tially even though the physics is entirely local. In such 
cases, however, it is expected that the spatial structure 
of Sxi^jt) will be trivial, and the condition Xo ^ ^ (in- 
dependent of temperature) will be violated. In computer 
simulations, these conditions may be checked, as we do 
in Sec.V. 

When the relaxation time diverges in a Vogcl-Fulcher 
manner, i.e. Tq = tq exp[DTo/(r — To)], one finds that 
the corresponding dynamic correlation volume also di- 
verges at To, as: 

DTTq 2 
VT(Taj ^ -77?, OC (IuTq) , 



(21) 



(T - ro)2 

where the last estimate holds sufficiently close to To. 

More generally, one can study the behavior of 
Xt(0, t) ~ 9Co(0, t; T)/dT as a function of time. Since at 
ah temperature Co(0, t = 0; T) = 1 and Co(0, oo; T) = 0, 
it is clear that XT(0,i) is zero at short and long times. 
We illustrate in Fig.[l]the shape of xt(0, t) for two glass- 
formers studied by molecular dynamics simulations de- 
scribed in Sec. El It has a peak for t « Tq. It is a useful 
exercise to study the example where the correlation func- 
tion is a stretched exponential with exponent (3 [not to 
be confused with l/fc^T], in which case: 



dCo{'Q,t-T) dlnr. 



91nT 



d In T ^ \Tn 



exp 



(22) 



This function behaves as a power-law, ^ at small times 
and reaches a maximum for t = Tc, before decaying to 
zero. The power-law at small times appears in the con- 
text of many different models, as discussed for the time 
behavior of X'i{i) Note also that for t = and 



T = Tg, one has: 



aa(0,T„;T) 



ainT 



l3m In 10, 



(23) 



where 



= Td\ogiQTa/dT\j, is the steepness index, 
which characterizes the fragility of the glass. Note that 
in many cases, the resulting numerical value of vt cc xt 
turns out to be already large in the late /3-regime, mean- 
ing that the concept of a cage is misleading because 
cagin g if fact involves the correlated motion of many par- 
ticles (H,!!^. 

Using the inequality in Eq. (fT^ with the results of the 
present section, we finally obtain a lower bound on the 
dynamical susceptibility X4(0 foi' Newtonian systems in 
the NPT ensemble, which is experimentally accessible: 



X4(r,t) > 



CP 



_1_ f dCo{r,t;T) 
CP V dhiT 



(24) 



This bound implies that as soon as xt increases faster 
than T~^ at low temperatures, X4 "^iH eventually exceed 
unity; since Xi is the space integral of a quantity bounded 
from above, this again means that the lengthscale over 
which the four-point correlation S4{y,Ta) extends has to 
grow as the system gets slower and slower. More quan- 
titative statements require information on the amplitude 
and shape of Si{y, To) which general field-theoretical and 
numerical results provide. 

The above result in Eq. ((M)) is extremely general and 
applies to different situations discussed in the next sec- 
tion. It however does not apply when the dynamics is 
not Newtonian, as for instance for Brownian particles or 
in Monte-Carlo numerical simulations [4l|, [1^, . The 
reason is that in these cases, not only the initial probabil- 
ity but also the transition probability from the initial to 
the final configuration itself explicitly depends on tem- 
perature. In Brownian dynamics, for example, the noise 
in the Langevin equation depends on temperature [55| . 
Hence, dCo{r, t; T)/dT receives extra contributions from 
the whole trajectory, that depend on the explicit choice 
of dynamics. We will argue below that when a dynami- 
cal critical point exists or is narrowly avoided, a system 
with Brownian dynamics should display dynamical cor- 
relations of the form X4 ~ Xt rather than the scaling 
Xi ~ Xt suggested by the above bound, Eq. 



D. Several generalizations 

1. Density rattier ttian temperature 

In the above section, we have shown that the response 
of the correlator to a change of temperature is related to 
dynamic correlations. Other perturbing fields may also 
be relevant, such as density, pressure, concentration of 
species in the case of mixtures, etc. For example, for 
hard-sphere colloids, temperature plays very little role 
whereas small changes of density can lead to enormous 
changes in relaxation times [STj . Using the expression 
for the probability of initial configurations in the NPT 
ensemble, and the fact that the dynamics only depends 
on the initial condition, one now derives the following 
equality: 



dCoir,t;P) 



dP 



ksT 



dx (o(x + r,i)o(x,0)(5w(0,0)), 

(25) 

which can again be used to define a dynamic correlation 
volume xp- Introducing the isothermal compressibility 
kt = {dp/dP)\T/ Pq and noting that the total variance 
of volume fluctuations per particle is given by ksTnT/ Poi 
we find: 



NEcv = Po 

= —ksTnT 



J dx (o(x + r, t)o(x, t = 0)Sv{0, t = 0)) 
dCo{r,t;p) 



d\np 



(26) 



7 




FIG. 1: (a) and (b) respectively show the self-intermediate scattering functions Fs(k,t) as a function of time for various 
temperatures in a binary Lennard- Jones mixture and the BKS model for silica, obtained from the molecular dynamics numerical 
simulations discussed in Sec.|V] (a) T = 2.1, 2.0, 1.05, 1.0, 0.75, 0.72, 0.61, 0.6, 0.51, 0.5, 0.47, 0.46, 0.435 and 0.43 from left 
to right, (b) T = 6100, 5900, 4700, 4600, 4000, 3920, 3580, 3520, 3250, 3200, 3000, 2960, 2750, and 2715 K from left to right. 
The arrows illustrate how xt ~ dFs/dT is obtained by finite difference for each pair of temperatures, (c) and (d) show the 
resulting xt(^) for both models, normalized by the strength of energy fluctuations. We show the absolute value since xt is a 
negative quantity. In both liquids the dynamic susceptibility presents a peak for t ~ Tc whose height increases as temperature 
decreases, revealing increasingly heterogeneous and spatially correlated dynamics. 



from which we deduce a second bound on the dynamic 
correlation volume X4(t): 



X4(r,*) > PaksTKT 



r dCo{r,t;p) 
\ 9 In p 



(27) 



Again, the right hand side of this expression is accessible 
to experiments [STj . Very importantly, and contrarily 
to the case of temperature, this inequality holds even 
for Brownian dynamics since the statistics of trajecto- 
ries have no explicit dependence on pressure or density. 
Finally, a similar inequality holds for binary mixtures, re- 
lating Xi {t) to the dependence of the correlation function 
on the mixture composition. 



2. Correlation and response in frequency space 

We have considered up to now the variance of the cor- 
relation function for a given time i, related to the four- 
point susceptibility Xiit)j but this can be generalized to 
the covariancc of the correlation in frequency space. 



Defining C'o(r,a;) = dtcos{ujt)Coir,t), the fluctu- 
ations of Co(j,uj) define a four-point susceptibility in 
Fourier space X4{r,uj) given by: 



cc 



(28) 



Repeating the same argument developed for correlations 
in time space, one finds: 



1 f dCo{r,Lo;T) 



(29) 



We have up to now considered correlation functions, but 
the very same string of arguments also applies to linear 
response functions, which can, in the context of New- 
tonian dynamics, be written solely as functions of the 
initial condition. For example, the susceptibility of the 
observable o to an external field X is: 



Xoir,t) 



1 



6o{x + r, t) 



exp 



-13 j dx' h{yi',t = 0) 



(30) 
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from which all the above results, transposed to response 
functions, can be derived. This is an important remark, 
since response functions, such as frequency dependent 
dielectric response or elastic moduli, arc routinely mea- 
sured in glassy materials. Their temperature or density 
dependence is therefore a direct probe of the dynamic 
correlation in these materials 1371. 



3. Higher derivatives 

One can of course study higher derivatives of the corre- 
lation functions with respect to temperature, which lead 
to higher order multi-point correlations between dynam- 
ics and energy or density fluctuations. For example, the 
second derivative gives a connected four-point correlation 
function: 



= J d^dy (o(x, t)o(x, 0)5e(y, 0),5e(0, 0))^. 

(31) 

The right hand side now defines a squared correlation 
volume, where the left hand side, computed for t = Ta, 
contains terms proportional to luTa/dlnT^ and to 
{d\nTa/d\nT)^. In most cases where Iutq diverges as 
an inverse power of temperature, or in a Vogcl-Fulchcr- 
like manner, one finds that the latter term dominates 
over the former. This means that this squared correla- 
tion volume in fact behaves like Xt- The same argument 
also holds for higher derivatives. 



E. Fluctuations and ensembles 



Because this result is important throughout this pa- 
per, we sketch here its proof, using ideas and a notation 
which should make clear the analogy with a similar result 
derived in Sec. IIIII using a field theoretical representation 
for the dynamics of supercooled liquids. Without loss of 
generality, we can choose the mean of all z^'s to be zero. 
The unconstrained joint distribution of the z^'s can be 
written as: 



VdetD 



(33) 



where _D is a certain M x M symmetric positive definite 
matrix. The unconstrained covariancc between Zi and Zj 
is well-known to be given by: 



(34) 



Let us now write D as blocks corresponding to the (m — 1) 
fluctuating variables and the (M — m -1-1) fixed variables: 



D = 



A B 

St C 



(35) 



where A is (m — 1) x (m — 1), _B is (m — 1) x (M— m + 1) 
and C is (A/ — m + 1) x (A/ — to+ 1). When the variables 
Zm, ...,ZAf are fixed, the unconstrained variables acquire 
non-zero average values which are easily found to be given 
by: 



M 



(36) 



Constrained vs. unconstrained fluctuations 



The above upper bounds in Eq. ((241) can in fact be 
given a much more precise meaning by realizing that fluc- 
tuations of thermodynamic quantities are Gaussian in the 
large volume limit [58| . except at a critical point. This 
allows one to show the following general result. Consider 
an observable O that depends on M Gaussian random 
variables zi, Z2, zm- We want to compare the ensem- 
ble where all the z^'s are free to fluctuate with the ensem- 
ble where one constrains a subset of the z^, say Zm, zm 
to take fixed values, with no fluctuations. In the limit of 
small fluctuations, the variances of O in the two ensem- 
bles are related through: 



M 



, Zm)c + 



a(o|z„,...,ZA/)a(o|z„,...,ZA/) , , , 

~ ~ {ZaZij)ciS2) 



a,f3—7n 



dz 



dz(3 



where the average in the ensemble where ZmT--,ZM are 
fixed is denoted by {■\zjn, Zm)- The subscript c means 
that we consider connected averages and we use Greek 
indices for the {M — m + \) constrained variables. 



To establish the relation between constrained and uncon- 
strained covariances, we note the following block matrix 
inversion rule D^^ =: 



{A- BC-^B^-^ -{A~ BC-^B^}-^BC-^ 
~C^^B^{A- BC-^B'^}-^ {C-B^A-^B}-^ 

(37) 

together with the matrix identity: 



{A- BC-^B^}-^ = A-^ 



(38) 



{A-^B){C - B^A-^By^{A-^B)^ 

The constrained covariancc (z^z^ |zm, zm)c is clearly 
given by {A~^)ij. Using the above identities, we directly 
obtain: 



(ZjZj) = {ZiZj\z^, Zm)c + '^^^p-^^{ZaZi3) 

^ aza CIZp 



(39) 



Now, the final result Eq. ([5^ above can be established 
simply by considering, to lowest order in the fluctuations, 
the observable O as an Af -|- 1th Gaussian variable cor- 
related with all the z^'s and apply the above equality to 
i = j = M + \. 
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2. From NPH to NPT 



4- Various sources of fluctuations 



Let us apply the general result Eq. ([321) to the ease of 
interest here, first to the case M = 1, with zi = H and 
number of partieles fixed. The two ensembles correspond 
to NPH and NPT, respectively. The above formula 
can be used with the correlation Co as an observable 
provided the dynamics is conservative, as argued above. 
Therefore: 



Xr^(r,i)=X4™(r,t) + ^(^^^^^') , (40) 



where we have replaced in the second term in the right 
hand side d/dH by {l/NcpkB)d/dT; xf^"ir,t) is the 
variance of the correlation function in the NPH ensem- 
ble where enthalpy does not fluctuate, a manifestly non- 
negative quantity. Therefore, the above equation recov- 
ers the lower bound Eq. ([M)) . with a physically explicit 
expression for the missing piece. The relative contribu- 
tion of the two terms determining Xa^'^ will be discussed 
in concrete cases in Sees. IIIll andlVl 



3. Local vs. global fluctuations 

The above discussion may appear puzzling for the fol- 
lowing reason: we have seen that the susceptibility x^it) 
is the space integral of a four-point correlation function 
S'4(y,t) which, although developing some spatial corre- 
lations on approaching the glass transition, remains rel- 
atively short-range in the supercooled liquid phase and 
should not depend on far away boundary conditions that 
ultimately decide whether energy is conserved or not. 
Since 6*4 (y,t) does not depend, in the thermodynamic 
limit, on the ensemble, how can its integral over space, 
X4{t), be affected by the choice of ensemble? The answer 
is that while the finite volume corrections to S'4(y,t) for 
a given y tend to zero when V 00, the integral over 
space of these corrections remain finite in that limit [s^ , 
and explain the difference between Xi^'^ 3'iid Xi^^- 
We understand that the physical correlation volume is 
given by Xa^'^'-, the long-range nature of the fixed en- 
ergy constraint leads to an underestimate of X4 in the 
NPH ensemble, which is irrelevant to describe local cor- 
relations. This is particularly important in numerical 
simulations [111: the study of S4{q,t) (the Fourier trans- 
form of 5*4 (y,t)) in the microcanonical ensemble will 
lead to a singular behavior associated to the fact that 
liniq^o 7^ 'S'4(q — 0,t), whereas the two coincide 

only in the ensemble where all conserved quantities are 
free to fluctuate {NPH for monoatomie liquids). The for- 
mer quantity is the physical quantity independent of the 
ensemble and will be denoted lim^^o 'S'4(q, t) = xX in the 
following, whereas the latter depends on the macroscopic 
constraint. We summarize this important discussion in 
SecHVl 



Equation ([32]) makes precise the intuition that dy- 
namic fluctuations are partly induced by the fluctuations 
of quantities that physically affect the dynamic behav- 
ior [U, Among these quantities, some are conserved 
thermodynamic quantities, such as the energy or density, 
and the dependence of the dynamics on those quantities 
are simply measured by the derivatives of the correlation 
function. The contribution of the local fluctuations of 
these quantities can therefore be estimated and lead to 
a lower bound to the total dynamic fluctuations. In a 
supercooled liquid one expects on general grounds that 
energy and density should play major roles in the dynam- 
ics. From the thermodynamic theory of fluctuations [60j . 
we know that in fact temperature (seen formally as a 
function of energy and density) and density are indepen- 
dent random variables, with variance {ST"^) = T"^ /{Ncy) 
and ((5u2) = ksTKr/iNpo). Therefore Eq. §^ gives for 
the "true" dynamic susceptibility: 



* 1 

X4 = — 
Cv 



dCo 
dlnT 



poksTKr 



dCo 
dlnp 



(41) 

The question of whether other, "hidden" variables also 
contribute to the dynamic fluctuations is tantamount to 
comparing X4^^ with xl- This question is very difficult 
to resolve theoretically in general. The rest of this paper 
and the companion paper [s^ are devoted to theoreti- 
cal arguments and numerical simulations which attempt 
to clarify this issue. Our numerical results suggest that 
^^VE ^ igg^g^ close to the glass transition, but 

that both X4^^ ^'^d X4 a-re in fact governed by the very 
same physical mechanism and define the same dynamical 
correlation length. 

Whether energy or density fluctuations is the domi- 
nant factor can be assessed by comparing the two ex- 
plicit terms appearing in the right hand side of Eq. ([4ip . 
Assuming time-temperature superposition, the ratio r of 
the two terms for t = r„ reads: 



= pocykBTKr 



d In Tn I ■ 
d\np 
d In Tq I 
dlnT IP , 



(42) 



Following Ref. [6IJ, and noting that poCvksTnT < 1 in 
usual liquids, we conclude that for most glass-formers, r 
is significantly less than one, which means that density 
effects are weaker than temperature effects and conse- 
quently contribute little to dynamic fluctuations. The 
situation is of course completely the opposite in hard- 
sphere colloidal glasses, where dhiTa/dliiT\p —>■ and 
r » 1. 



F. Summary 

After motivating the use of multi-point correlation 
functions to detect non-trivial dynamic correlations in 
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amorphous materials, we discussed the idea that induced 
fluctuations are more easily accessible experimentally 
than spontaneous ones, and can be related to one an- 
other by fluctuation-dissipation theorems. Elaborating 
on this idea, we have shown that the derivative of the 
correlation function with respect to temperature or den- 
sity directly gives access to the volume integral of the 
correlation between local energy (or density fluctuations) 
and dynamics. This relation can be used to show on very 
general grounds that a sufficiently abrupt slowing down 
of the dynamics must be accompanied by the growth 
of a correlation volume. The detailed relation between 
these susceptibilities and a correlation lengthscale how- 
ever depends on the amplitude and spatial structure of 
the multi-point correlation functions. 

We have then shown that the dynamic four-point sus- 
ceptibility at q ~ Q, which corresponds to the fluctuation 
of global intensive dynamical correlators, depends in gen- 
eral on the chosen statistical ensemble. In the case where 
conserved variables are allowed to fluctuate, we showed 
that the dynamic four-point susceptibility is bounded 
from below by terms that capture the contribution of en- 
ergy and density fluctuations to dynamic heterogeneities. 
Our central results, suggesting a way to estimate a dy- 
namic correlation volume from experiments, are given in 
Eqs. ((Ml BT|) . Whereas we expect that for most su- 
percooled liquids, the contribution of temperature is the 
dominant effect, the quality of our bounds as quantitative 
estimators of X4: and their physical relevance is, at this 
stage of the discussion, an open question which we care- 
fully address below, in particular in Sec. IIVI and in the 
companion paper [38[ . The following section is devoted to 
a quantitative study of this question within a field-theory 
formalism. A surprising outcome of this analysis is that 
the dynamic four-point susceptibility at g = correla- 
tions not only depend on the chosen statistical ensemble, 
as shown above, but also on the choice of microscopic dy- 
namics whether Newtonian or stochastic. Of course the 
dynamic four-point susceptibility at non-zero q depends 
only on the choice of microscopic dynamics. 



III. CORRELATION OF DYNAMICAL 
FLUCTUATIONS: A FIELD-THEORETICAL 
PERSPECTIVE 



thermore predicts that, contrary to the behavior of cor- 
relators measuring the average dynamics, the details of 
dynamic fluctuations depend on the dynamics in a re- 
markable way. However, since in all cases, the object 
responsible for the increase of these dynamic correlations 
is the very same susceptibility A~^, the physics revealed 
by the correlations is independent both of the ensemble 
and of the dynamics, and genuinely reflects the collective 
nature of glassy dynamics. 

In the companion paper [sst we will point out how sim- 
plifications can occur if a true dynamical critical point 
exists, as within mode-coupling theory, a particular self- 
consistent resummation scheme. In the following we aim 
instead at keeping the discussion more general than the 
confines of mode-coupling theory or any other particu- 
lar theoretical approach. This is important since mode- 
coupling theory is not expected to apply close to the glass 
transition temperature, whereas the present physical con- 
clusions do. 



A. The dynamic field-theory 

1. A reminder of the usual static case 

The dynamic field-theory strategy is analogous to the 
one used for ordinary static critical phenomena which we 
now recall, focusing on the ferromagnetic Ising transition 
as a pedagogical example [g^]. The starting point is the 
Legendre functional transform r(m(x)) of the free energy 
/3F(/i(x)), itself defined as a functional of the magnetic 
field /i(x): 



r(m(x)) = PF{h{x)) - / dx'/i(x')m(x'), 



(43) 



where h{x.) on the right hand side is the field that leads 
to the magnetization profile 771 (x). The magnetization is 
determined via the equation: 



m(x) 



Sh{x) ■ 



(44) 



Two important properties of the functional r(m(x)) that 
can be directly derived using the previous relation are: 



In the following we develop in detail an approach to 
dynamical fluctuations in supercooled liquids based on 
general field-theory techniques, and discuss how a non- 
trivial lengthscale can be generated by interactions and 
manifests itself in quantities like Xi or XT- Wc identify 
precisely the 'susceptibility' (called below) respon- 
sible for all interesting dynamic correlations. We dis- 
cuss the origin of the ensemble dependence of dynamic 
fiuctuations described above from a diagrammatic point 
of view. This is important since any self-consistent re- 
summation or approximation scheme must be compati- 
ble with the bounds derived above. This formalism fur- 



ST 



(5m(x) 



(.(x).(x'))]-^ (45) 



The last exact identity indicates that the operator ob- 
tained by differentiating the functional T twice is the in- 
verse of the spin-spin correlation function (considered as 
an operator). Note that these are simple generalizations 
of usual thermodynamic relations. 

In general one cannot compute F exactly, but one can 
guess its form using symmetry arguments, and compute 
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it approximately in a perturbative (diagrammatic) ex- 
pansion in some parameter. Using the above identities, 
no further approximation is needed to obtain correla- 
tion functions. In its simplest version, F corresponds to 
the Ginzburg-Landau free energy functional. The sad- 
dle point equation for the magnetization then leads to 
the mean-field description of the transition, whereas the 
second derivative term gives the mean-field result for the 
spin-spin correlation function, valid when the space di- 
mensionality is sufficiently large. 

In the following, we will present a theory of dynamic 
fluctuations within a field-theoretic framework similar to 
the above static formalism. The main difference is that 
in the context of glassy dynamics, the relevant order pa- 
rameter is no longer a one-point function like the magne- 
tization but instead a two-point dynamic function which 
has to be introduced as an effective degree of freedom in 
the dynamic free-energy functional. 

2. Dynamic free-energy Junctionals and fluctuations 

Different dynamic field-theories have been used in the 
literature to analyze the dynamics of dense liquids. The 
common strategy is to write down exact or phenomeno- 
logical stochastic equations for the evolution of the slow 
conserved degrees of freedom. For instance, for Brow- 



where S is the action of the field theory, ha 's are such that 
($a) = and Ka,b imposes a certain value for the two- 
point functions: {^a^b} — ^a'^b = Ga,b- The properties 
of r(\E'Q, Ga,b) are the same as in the static case because 
formally it is the same mathematical object. The only 
difference is that the dynamical functional depends on a 
larger number of variables. The difficulty is to devise an 
approximate expression for the functional F. Once this 
is done, one should differentiate the functional once to 
obtain self-consistent equations for the order parameters 
Ga,b and twice to obtain (after inversion) and expres- 
sion for their fluctuations. More precisely, we introduce 
the following matrix of second derivatives: 



ST 


ST 1 




■ A 


B ' 


SGa.bSGa.d 


SGa b<5*e 




ST 


ST 




St 


C 


S-HjSGa b 







where we have introduced three block matrices A, B, C, 
in full correspondence with those introduced above in 



nian dynamics the only conserved quantity is the local 
density (energy and momentum are not conserved). The 
equation for the local density is in that case the so- 
called Dean-Kawasaki equation [6^, [g^l, which can be 
derived exactly for Langevin particles (see Refs. [1^ [6^ 
for a discussion of different field-theories associated with 
such dynamics). In general, the field-theory associated 
to a given stochastic dynamics is obtained through the 
Martin-Siggia-Rose-deDominicis-Janssen method, where 
one first introduces response fields enforcing the cor- 
rect time evolution and then averages over the stochastic 
noise 

We will use a general notation that will allow us to 
treat all field theories proposed in the literature [gI, 
on the same footing. In all those field theories one has 
a set of slow conserved fields, (jii {i = 1,...,to), and 
the corresponding response fields, ipi arising from the 
Martin-Siggia-Rose procedure [b^]- It will also be use- 
ful to put 4>i, (pi into a single 2m dimensional vector $a, 
a = l,...,2m. The average over the dynamic action of 
<^a will be denoted ^a' (^a) = ^a- As in the static case, 
the starting point of the analysis is a Legendre functional 
(also called the generator of two-particle irreducible di- 
agrams or Baym-Kadanoff functional) (g^, [H, [6^ . It is 
equal to: 



(46) 



(47) 



Sec. Ill E II The inversion of d^V allows one to obtain 
the objects of interest in this paper. For example, inver- 
sion in the "GG-sector" defines the four-point space-time 
correlation functions: 

{d'TVaMld = (($a(x,t)<i>,(x',t') - Vl/,(x,t)*,(x',t')) 

X (<i>c(y, s)$<i(y', s') - *,(y, s)-^a{y\ s')))c, (48) 

where (•)c means that we are focusing on the connected 
component. Similarly, inversion in the "G^'-sector" 
defines the three-point functions, such as the energy- 
correlation correlator defined in the previous section, 
whereas inversion in the "^"I'-sector" leads to the ex- 
act propagators of the conserved quantity. For example, 
when 5" is the energy, one obtains the exact energy prop- 
agator (dressed by interactions) , which is expected to be 
diffusive in the hydrodynamic limit. 

At this stage, it is important to recall that the dynam- 
ical functional F has a direct diagrammatic expression 



F(*a,Ga,b) = -In j V^a CXp {^S{{<^a})- j dtdx /l^ (x, <) (x, t) - *a (x, t)] 

-- / dtdt'd^d^i' Ji:a,6(x,f;x',i')[*a(x,i)$6(x',t')-*a(x,i)*b(x',t')-Ga,6(x,i;x',i') 



a,b=l 
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r(*,G) = --TrlogG+-TrGoi[G + **]-$2P/(*,G), 

(49) 

where $2P/(^,G) is the sum of aU two particle irre- 
ducible Feynman diagrams (that cannot be decomposed 
in two disjoint pieces by cutting two lines) constructed 
with the vertices of the theory and usin g th e full prop- 
agator G as lines and as sources [6^ l68l. [69| . Both 
the internal indices and spatio-temporal arguments were 
skipped for simplicity. The first derivatives lead to the 
self-consistent equations for the order parameter. Since 
in dynamical field-theories for liquids the slow physi- 
cal fields are in fact conserved quantities the equations 
5T/5'^a = do not fix the values of the physical fields 
that have to be fixed by the initial conditions. On the 
other hand; they set to zero the average of the response 
fields and enforce translational invariance (soj . 

On the other hand, the derivatives ST/SG = lead to 
formally exact self-consistent equations for the two-point 
correlation functions. These equations can be written as 
a Schwinger-Dyson matrix equation: 

G-i=Go-i-S(G), l](G) = ^, 

where S is the self-energy. A given approximation con- 
sists in retaining a given set of diagrams in $2P/, or al- 
ternatively in S(G). For example, mode-coupling the- 
ories generically consist in only retaining the "bubble" 

diagram for S(G) , see Refs. m m, [zs 

for detailed dis- 
cussions, and [38| in the present context. 

B. Three-point correlation: dynamic susceptibility 
and hydrodynamic contributions 



n 



z 

n 




FIG. 2: Diagrammatic representation of the parquet dia- 
grams obtained by expanding Eq. (|53p : (1 — dc'^GG)'^ = 
E„(9gEGG)". 

read 

5t = x*(**>, (52) 

providing an exact decomposition with a simple physical 
meaning. The correlation between energy at one point 
in space-time and dynamics elsewhere is governed by the 
sensitivity of the dynamics to energy changes, as encapsu- 
lated by X*, which contains all genuine collective effects 
in the dynamics induced by interactions. This correla- 
tion is mediated by energy transport, (^'5'), with has a 
trivial hydrodynamic structure. 

In order to see this more clearly, let us now explore 
the diagrammatic content of x* ~ A^^B. The three- leg 
vertex contribution B = sqs^ is generically expected to 
be non singular. The A^^ term, on the other hand, can 
be rewritten using the general expression of F as: 



From the above general inversion formulas for block 
matrices, Eq. (j37[) . one can obtain an expression of the 
inverse of d^T in the "G\l/-sector" in a form transparent 
both from physical and diagrammatic standpoints. The 
off-diagonal block-element of Eq. ([57)1 gives, in particular, 
the energy-dynamics correlator (see Eq. ^6]) ) and can be 
rewritten exactly as: 

= (9'r)-i'^* = (50) 

where we have used that (^"5) = (O^F)"!^**. Now, 
the equation determining the two-point correlators is 
6r/6G = 0. Therefore the variation of the value of G 
due to a small variation of ^, all other parameters being 
kept fixed, is given by: 



(5* 



5^r 



SGSG 



SGS'^i'^ 



(5* 



-A-^BS"^, (51) 



showing that the operator = —A^^B is the response 
of two-point correlators to a change in conserved quan- 
tities. Gathering these results, three-point functions 



A- 





-1 


SGSG 





where the objects in the above expression are four-index 
matrices. This term can be rearranged as follows: 



SGSG 



,,b:c,d 



Ga,c'Gb,d' 



^ I ^^^^ 

c" .d" ' 



One can now formally expand the term in parenthesis 
as (1 - 9gSGG)-i = E„(^SGG)" to recover the so- 
called "parquet" diagrams [631 ^^^^^ S^'^^ a formally exact 
representation of the four-point function, see Fig.O This 
infinite series can provide a divergent contribution (as is 
the case within MCT at the critical point [l3|)j signaling 
the existence of a growing dynamical correlation length 
and non-trivial collective effects. 
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The important conclusion of this section is that the 
three-point function contains both a long-ranged hydro- 
dynamical contribution (5'\E') related to energy conser- 
vation, and an interaction specific contribution - the dy- 
namical susceptibility . When specializing to the inte- 
gral over space of the three point function, as in Eq. (fT6|) . 
the contribution of (4'^') factors out and gives thermody- 
namic prefactors. This shows that dC /dT gives in fact a 
direct access to the dynamical susceptibility x.f at 9 = 0. 
Therefore, the lengthscale extracted from dC /dT reveals 
the existence of collective dynamics, and is not related to 
any thermal diffusion or other hydrodynamical length. 



C. Four-point correlation functions and ensemble 
dependence 

Let us now turn to a similar analysis of the four-point 
correlations. We start again from the above general in- 
version formulas for block matrices, Eqs. ([571155)1 . In the 
simple case where the '^aS are identically zero by sym- 
metry, as happens for instance in the p-spin model for 
which a gauge symmetry implies that the average value 
of the spins is always zero, the block matrix B is also 
zero. Equation ([37]) then simplifies to: 



n -1 



SGa,bSGc,d 



A' 



(54) 



In general this symmetry does not hold, in particular for 
liquids for which the analysis is more involved. However, 
it turns out that remains the fundamental object. 
Using Eq. ([55)1 and the bottom right part of the ma- 
trix inversion relation, Eq. ()37p . the four-point correla- 
tion functions can be written in a physically transparent 
way: 



(S'r)-!;^, 



6^r 



SGa,bSGc,d 
SGab 



f6G. 



cd 



(55) 



This expression parallels Eq. ([4T)l in Sec. |TT1 and the last 
term corresponds to the dynamic fluctuations induced by 
the fluctuations of conserved quantities. This formula is 
however much more general because it applies not only 
to XiW but also to 54 (g, t). Indeed, in Fourier space, the 
terms contributing to Sii{q,t) read: 

Therefore the extra contribution from conserved quanti- 
ties, namely the last term in Eq. (|55p . reads: 



E 



d{5p.uM^pu,+m 

9*e(c^,g) 



:«'e(w, (?)*/( 



r 



-1)) 



d{Sp^kAt)Spk,-gm 



(57) 



Now, one should notice that all terms corresponding to 
indices in the response field sector of ^' (i.e. e, / > m) 
identically vanish at g = 0. The reason is that the re- 
sponse fields always appear in the vertices of the field- 
theory in the form V^*. As a consequence, terms like 
S^r /6G6'$e{^, q) , for e > m, are proportional to q at 
small q. 

In the case q = 0, the value of conserved fields such as 
\['e(a;,q') for e < m are by definition constant over time 
and set by initial conditions: (^e(tj, (7)^'y (— cj, — (j))c = 
VS{uj)Tief where Se/ are the correlators of thermody- 
namic fluctuations of all conserved quantities 'J, deter- 
mined by the probability distribution of initial condi- 
tions. As a consequence the term in Eq. ([FT)) at g = 
precisely reduces to the form discussed in the previous 
section on general grounds for Xii't) = '5'4('? = 0, t): 

^ ■ ^ ' 



9* 



/ 



For Brownian dynamics, density is the only conserved 
quantity and thus only one term, m = 1, contributes to 



the sum in Eq. ([58p . In the case of Newtonian dynamics 
there are in principle 2 + d conserved quantities, density, 
momentum and energy. However, by symmetry, the con- 
tribution of the momentum fluctuations is zero, and only 
density and energy should be considered. In p-spin dis- 
ordered systems, on the other hand, this extra term is 
absent, and m = 0. 

The conclusion is once again that the choice of statis- 
tical ensemble matters for determining fluctuations of in- 
tensive dynamical correlators which correspond to q = 0. 
For q 0, the extra terms Eq. ()57p are in general al- 
ways non-zero and contribute to S^q^t). On the other 
hand, if one focuses on the case where q = exactly, 
the initial distribution is crucial. As an example, in the 
case of Newtonian dynamics in the NVE ensemble, all 
the extra contributions vanish since in that ensemble all 
conserved quantities are strictly fixed and T,^^^ = 0. 
Thus, we find again within this formalism that the equal- 
ity limg_»o Si{q,t) = 84(0^ t) is valid only in the ensem- 
ble where all conserved quantities fiuctuates. In other 
ensembles, such as NVE, the limit is singular. 
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D. A direct measure of dynamical susceptibility 



n 
m 




5r _ 




5G5\|/ 



FIG. 3: "Squared-parquets" representation of the contribu- 
tion of conserved quantity fluctuations to X4- They corre- 
spond to the second term in Eq. (|55|) . 



Let us now explore the diagrammatic content of 
Eq. ([55)1 . The first term was already discussed in the 
previous section and can be expressed as a sum of par- 
quet diagrams, see Fig. O The second term in Eq. fSS]) 
also has a direct diagrammatic interpretation, shown in 
Fig. [31 It consists of two parquets closed by three- leg ver- 
tices and joined by a correlation function of conserved 
variables. The wavevcctor q in Si{q,t) is, in the dia- 
grams, the wavevcctor flowing into the parquets and in 
the middle "link" corresponding to (^'e5'/)c, as it ap- 
pears explicitly in Eq. ([57[) . A detailed analysis of the 
structure of the diagrams shows that hydrodynamic scal- 
ing between time and length is present only in the middle 
"link" corresponding to (\I'e^/)c. 

As discussed in [13], MCT provides a simple approx- 
imation in which the self-energy E is approximated by 
the "bubble" diagram (see (S^). Within this approx- 
imation the parquet diagrams simplify into the ladder 
diagrams analyzed in 17[, which diverge at the mode- 
coupling critical point. In Ref. fv^ . however, only the 
ladder diagrams were analyzed. The contribution to 
54 of Eq. ([57[) corresponding to the "squared ladders" 
was overlooked. As a consequence the MCT results in 
Ref. [l3| for X4{t) = ^^{q = 0,t) only apply close to 
the critical point in the following cases (see also [s^ for 
further discussion): 

• NVE ensemble for Newtonian dynamics; 

• NVT ensemble for Brownian dynamics; 



• p-spin models. 

On the other hand, whenever conserved quantities arc 
allowed to fluctuate, or when considering Si{q^ t) at non- 
zero values of g, the contribution of ([57p may be impor- 
tant. For example, within the context of MCT where 
diverges as (where e is the reduced distance from the 
critical point), the contribution of ([55)) in fact becomes 
dominant and Xi foi' Newtonian dynamics diverges much 
faster, as e~^. However, see II for a discussion of the 
application of these MCT results to real systems, where 
the MCT transition is avoided. 



The analysis of the above sections show that in gen- 
eral Si{q,t) and X4(0 receive contributions of different 
physical origin with possibly different temperature de- 
pendencies, and whose relative amplitude might even de- 
pend on the chosen microscopic dynamics (Brownian or 
Newtonian). On the other hand, we have seen that all 
the interesting physics is contained in the fundamental 
operator A = / 5G5G], which governs the growth of 
dynamic correlations. Therefore it is both of theoretical 
and practical importance to introduce an observable with 
a physical content similar to that of S4{q^t), but unaf- 
fected by the presence of global conservation laws and 
therefore by the choice of statistical ensemble. Such an 
observable was discussed recently HI] . It corresponds to 
the response of the intermediate scattering function (the 
two-point correlator G) to a small inhomogeneous exter- 
nal potential Vext- Within the previous formalism one 
writes: 



S^T SG 
SGSGSvZf 



6G6V,^ 



SGS^WZt 



= 0, 



and therefore: 



SG 

SVexi 



\5G5Vext SGS'^SVe, 



Since the source term on which the operator acts 
is expected to be only weakly temperature and density 
dependent, one sees that this quantity gives an almost di- 
rect measure of the critical behavior of the dynamic cor- 
relations encoded in the operator A. When the external 
potential is homogeneous in space one finds a quantity 
proportional to x* above [39| . while for an inhomoge- 
neous external potential one can probe the full spatial 
structure of dynamic fluctuations. Indeed, when one dif- 
ferentiates with respect to the Fourier component Vextiq)i 
the wavevector q plays the same role as for S4 [39j . This 
can be seen at the diagrammatic level because q is the 
wavevector entering into the ladders in Fig. [31 



IV. PHYSICAL CONSEQUENCES AND ISSUES 

At this stage, it is important to summarize the con- 
clusions drawn from the rather dense theoretical analysis 
presented above. This will allow us to identify clearly the 
questions that need to be tested numerically before possi- 
bly extrapolating these conclusions to real glass-forming 
systems. 

We established in the previous section that all non- 
trivial collective dynamical effects are encoded into a cer- 
tain operator A^^, which could in principle be reached 
by measuring the sensitivity of the local dynamics to an 
external potential [s^. More easily accessible quantities 
are derivatives of two-point correlations with respect to 
temperature, or density. We have shown in detail how 
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these are indeed proportional to and provide lower 
bounds on X4, and are therefore of direct interest to probe 
the growth of a dynamic lengthscale in glasses, as claimed 
in [37| . However, the assumption that growing suscepti- 
bilities imply growing lengthscales needs to be discussed 
more thoroughly. 

A. Growing susceptibilities vs growing lengthscales 

The first important remark is that the lower bound 
on X4 obtained in the previous sections is useful only 
when X4 is significantly larger than one because Xi is 
of the order one even in an ideal gas [U- The second 
remark is that one has to be sure that the growth of 
the susceptibility is due to a growing length and not to 
growing local fluctuations. For simplicity, suppose that 
the energy-dynamics correlator ST{y,Ta) can be written 
as S't(O) X Ks space integral xt is then 

given by: 

XT = STiO)e-" J dV-'-'^f{u). (59) 

This shows the origin of an increase in xt = dC/dT as 
T decreases is due to 

• either an increase in ^ with a roughly constant 
St{0) 

• or because S't(O) increases whereas ^ is trivial, 

or, of course, through a combination of both. In order 
to be confident that the first scenario is the correct one, 
and that xt can be used to estimate a correlation vol- 
ume, one needs to be sure that S't(O) is of order one 
and basically temperature independent. This requires in 
principle some extra information, for example on the full 
spatial dependence of S't(2/, Tq,). This will be checked in 
numerical simulations below. We note also that MCT 
precisely realizes the first scenario above. 

From a physical point of view, one expects the en- 
thalpy fluctuations 5h to contain a fast (kinetic) part 
and a slow (configurational) part, of similar order of 
magnitude (ksT). While it is clear that the fast part 
should have very small correlations with the local corre- 
lation on time scale Tq, there is no reason to think that 
{o{x,Ta)o{x,t = 0)Shsiow{x,t = 0)) is particularly small. 
Quite the contrary, we expect that this is of order ksT 
in glassy systems. But interestingly, this suggests that 
the specific heat cp that should enter the relation be- 
tween XT and ^ should be the so-called excess specific 
heat Acp, restricted to slow (glassy) degrees of freedom, 
as surmised in (s^ . 

B. Statistical ensemble and dynamics dependence 
of dynamic fluctuations 

A rather bizarre conclusion of the previous section is 
that global four-point correlators, corresponding to the 



fluctuations of intensive dynamical correlators, not only 
depend on the statistical ensemble (for g = 0) but, re- 
markably and perhaps unexpectedly, also on the choice 
of dynamics. This is to be contrasted with the case 
for two-point correlators, which are independent of the 
chosen ensemble and are known numerically to be inde- 
pendent of the dynamics, at least in the relevant "slow" 
regime [i^ |4l|, [7l|. This shows that four-point corre- 
lators, although containing some useful information on 
dynamical heterogeneities, mixes it with other, less inter- 
esting physical effects. Clear-cut statements with four- 
point quantities can however be made when the dynamic 
lengthscale grows substantially at some finite tempera- 
ture or density, as for example within MCT where the 
operator develops a zero mode that leads to a diver- 
gence of the dynamic lengthscale ^. When the dynamic 
lengthscale becomes very large, these statements may be 
summarized as follows : 

• X4(''a) for NVT Newtonian dynamics diverges 
more strongly than X4 ('''«) for NVT stochastic dy- 
namics; 

• X'ii'^a) for NVE Newtonian dynamics diverges like 
Xi{Ta) for NVT stochastic dynamics; 

• Xii'^a) for NVE Newtonian dynamics and NVT 
stochastic dynamics diverge like XT{Ta) (or 

Xp(^a)). 

We will test these statements numerically in the next 
section, and will indeed establish that xt, Xa^^ and xf 
increase in exactly the same way with Tq,. The full time 
dependence of these different correlators will be discussed 
in the companion paper, Ref. [38| . 



C. A unique dynamic correlation length 

Let us emphasize again that although X4(''q) for NVT 
Newtonian dynamics and stochastic dynamics diverge 
differently, our results strongly suggests these quanti- 
ties in fact reflect the same underlying physics, which 
is the growth of a unique lengthscale ^ in all of these 
cases. Only the relation between X4('''a) or XT{Ta) and 
^ changes: dynamic fluctuations are amplified because of 
conserved variables. This becomes clear when one con- 
siders the (ensemble- independent) function Si{q^Ta) for 
q ^ Q. In all of these cases, 6*4 (5, Tq) can be written 
as a scaling function 34 ((7^) with the same ^ but different 
functional forms. For example, (qS,) for Newtonian dy- 
namics can be written as gf{qC) + Cq\9i ^ where Cg 
is a coefficient and gf{q£,) ^ 9t{<iC} scaling function 
for Brownian dynamics or governing Sriq, To). Note that 
the relation between X4(''a) for NVE Newtonian dynam- 
ics, X4('''a) for NVT Brownian dynamics and XT{Ta) niay 
not be accurate far from any critical point, since these 
quantities are affected by different, non-critical prefac- 
tors. 
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D. Response vs correlation functions 



tides of types A and B, with interactions 



Four point correlators were originally hoped to be 
suited to quantify precisely dynamical heterogeneities in 
glass-formers, as motivated in Sec. [TTl The conclusion of 
the previous section and the numerical results of the fol- 
lowing ones show that although they contain indeed cru- 
cial information, it is mixed up with less interesting phys- 
ical effects. Nevertheless, a unique dynamic correlation 
length seems to govern the slowing down independently 
of the dynamics, global dynamic fluctuations depend on 
the dynamics and on the ensemble. As discussed formally 
in the previous section, response functions measuring the 
response of the dynamics to local perturbations do not 
present these difficulties. They should be independent of 
the microscopic dynamics, as it is the case for two-point 
correlators, and probe directly the dynamic correlations 
without mixing them up with other effects due to con- 
servation laws. 

In the following, we give some numerical evidence for 
the most important claims made in this paper: the exis- 
tence of unique lengthscale ^ governing the growth of Xi 
and XT and the ensemble and dynamics dependence of 
the four-point correlators. 



NUMERICAL RESULTS FOR TWO 
MOLECULAR GLASS-FORMERS 



We now present our numerical calculations of the dy- 
namic susceptibility XT{t), its relation to X4(0i ^^'^ 
the behavior of spatial correlations St and 6*4 in two 
well-studied models of molecular glass-formers: a bi- 
nary Lennard- Jones (LJ) mixture [23, considered as a 
simple model system for fragile supercooled liquids 
and the Beest, Kramer, and van Santen (BKS) model, 
which is a simple description of the strong glass-former 
silica [7^ [t^I. a first motivation for these simulations 
is that all terms contributing to the dynamic fluctu- 
ations can be separately evaluated and quantitatively 
compared. Spatial correlators and dynamic lengthscales 
can be directly evaluated in the simulations to conflrm 
the link between dynamic susceptibilities and dynamical 
lengthscales. Therefore, the claim made in Ref. [s^l that 
XT{t) yields direct experimental access to a dynamical 
lengthscale can be quantitatively established. A second 
interesting feature is that the influence of the microscopic 
dynamics and statistical ensemble can be quantified in 
the simulations by keeping the pair potential unchanged 
but by switching from the energy conserving Newtonian 
dynamics to some stochastic dynamics which locally sup- 
plies energy to the particles. 



A. Models and technical details 

The binary LJ system simulated in this work is a 80:20 
mixture of Na = 800 and Nb = 200 Lennard-Jones par- 



(- 
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r 



(60) 



where a, /3 G [A, B] and r is the distance between the 
particles of type a and (3. Interaction parameters eap 
and aaf3 are chosen to prevent crystallization and can 
be found in Ref. [t^I. The length, energy and time 
units are the standard Lennard-Jones units aAA (par- 
ticle diameter), saa (interaction energy), and tq = 

myicr^^/(48e^yi), where niA = ms is the particle mass 
and the subscript A refers to the majority species. Equi- 
librium properties of the system have been fully char- 
acterized [72|. At the reduced density po = 1-2, where 
all our simulations are carried out, the MCT transition 
has been conjectured to be in the vicinity of Tc « 0.435 
[t^ I. The slowing down of the dynamics, T > 0.47, can 
be correctly described by mode-coupling theory, but this 
description eventually breaks down when lowering the 
temperature further, T < 0.47 [t^ . 

To check the generality of our results, we have also 
investigated the behavior of a second glass-former, char- 
acterized by a very different fragility. To this end, we 
simulate a material with an Arrhenius dependence of its 
relaxation time, namely silica. Various simulations have 
shown that a reliable pair potential to simulate silica is 
the one proposed by BKS [tI, [tII . The functional form 
of the BKS potential is 



^ V Aai3eyiY>{~Bapr) f-, (61) 



where a, (3 G [Si, O] and r is the distance between the 
ions of type a and (3. The values of the constants 
qa,qi3, Aafs, Bap, and Cap can be found in Ref. [t^. For 
the sake of computational efficiency the short range part 
of the potential was truncated and shifted at 5.5 A. This 
truncation also has the benefit of improving the agree- 
ment between simulation and experiment with regard to 
the density of the amorphous glass at low temperatures. 
The system investigated has Ngj = 336 and Nq = 672 
ions in a cubic box with fixed size L = 24.23 A. The 
Coulombic part of the potential has been evaluated by 
means of the Ewald sum using a constant aL=10.177. 

For both LJ and BKS models we have numerically in- 
tegrated Newton's equations of motion using the velocity 
Verlet algorithm [sHi using a time-step h^'^ = O.OItq and 
j-iPKS _ 2.6 fs, respectively. Doing so we can measure 
spontaneous dynamic fluctuations in the microcanonical 
NVE ensemble. Before these microcanonical production 
runs, all systems are equilibrated using a stochastic heat 
bath for a duration significantly longer than the typical 
relaxation time, Tq,, implying that particles move over 
several times their own diameter during equilibration. 
Production runs were at least larger than 30to,, and sta- 
tistical convergence for dynamic fluctuations was further 
improved by simulating 10 independent samples of each 
system at each temperature. Repeating this strategy for 
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many temperatures in two molecular systems obviously 
represents a substantial numerical effort. 

To check the influence of the microscopic dynamics, 
and in particular the role of the energy conservation, 
we have also performed stochastic simulations of the 
LJ s yste m using two different techniques. Following 
Ref. [40[, we have simulated Brownian dynamics where 
Newton's equations are supplemented by a random force 
and a viscous friction whose amplitudes are related by 
the fluctuation-dissipation theorem. The numerical al- 
gorithm used to integrate these Brownian equations of 
motion is described in Refs. [iO, UHl using the time step 
of ft-BQ = 0.016ro and a friction coefficient C, = lOmro. 
We use the equilibrium configurations obtained by MD 
simulations as starting point for our production runs in 
Brownian simulations. Finally we have implemented a 
second stochastic dynamics, a standard Monte Carlo dy- 
namics, with the LJ potential [4l|. At time t, the particle 
i, located at the position ri{t), is chosen at random. The 
energy cost AE to move it to the new position (i) + S is 
evaluated, d being a random vector comprised in a square 
of lateral size Sj^ax = 0.15. The Metropolis acceptance 
rate, p ~ min(l, e"^^"^), is then used to decide whether 
the move is accepted [5a]. One Monte Carlo time-step 
represents N — Na + Nb attempts to make such a move. 

For BKS, we only present results for ND because BD 
simulations at low enough temperature would be numer- 
ically too costly in this system. The reason is that a 
very large friction coefficient is needed to have a truly 
damped dynamics [t^, making the overall relaxation 
much too slow to be studied numerically at low tem- 
perature. Monte Carlo simulations are similarly slow be- 
cause of the long-range character of the Coulomb term 
in the BKS potential. Very recently we have developed a 
short-range approximation of the BKS potential that al- 
lows much faster Monte-Carlo simulations, and we shall 
mention some preliminary results obtained for dynamic 
susceptibilities using this method [t^. 

B. Physical observables 

Following previous work p^ . [ssl . [3^ , we monitor the 
dynamical behavior of the molecular liquids through the 
self-intermediate scattering function, 

F,(k,t) = ^^£e'''-["-.W-r.(o)]^^ (62) 

where the sum in Eq. (|62p runs over one of the species 
of the considered liquid [A or B in the LJ, Si or O 
for silica). We denote by /s(k, t) the real part of the 
instantaneous value of this quantity, so that we have 
F,(k,t) = (/,(k,t)). 

The four-point susceptibility, X4(t), quantifies the 
strength of the spontaneous fluctuations around the av- 
erage dynamics by the variance, 

X4(i)=iVa[(/.'(k,i))-i^i^(k,t)] . (63) 



In principle, Xi{^) Eq. (j63p retains a dependence on 
the scattering vector k. Since the system is isotropic, 
we circularly average ([5^ and over wavevectors of 
fixed modulus. Note that the value of dynamical corre- 
lations depend on |k| as shown in [13, [s^. A detailed 
analysis of this dependence has been performed in (ssj 
and will be further discussed in [s^. In the following we 
will focus only on the value of |k| for which the dynam- 
ical correlations are more pronounced and that measure 
the correlation of the local dynamics. For the LJ sys- 
tem we will mainly consider results for |k| = 7.21 and for 
the BKS one |k| = 1.7 A~^. These values respectively 
represent the typical distance between A particles, and 
the size of the Si04 tetrahedra. As discussed above, we 
expect Xi{i) to depend on the chosen statistical ensem- 
ble, e.g. NVE or NVT, for Newtonian dynamics, and 
to depend also on which microscopic dynamics is chosen, 
stochastic or energy conserving. 

To evaluate the temperature derivatives involved in 



xHt) = —F,{k,t), (64) 



we perform simulations at nearby temperatures, T and 
T + ST, and estimate XT{t) through finite differencing, 
Xrit) ~ SFs{k,t)/ST, as illustrated by arrows in Fig. [1] 
in Sec. |TT1 For this procedure to be effective, temper- 
ature differences must be small enough that linear re- 
sponse holds. Taking ST too small leads however to poor 
statistics. The smallest ST which might be used can be 
estimated by comparing the statistical noise of Fs{k,t) 
to the expected response xt x ST. This leads in our case 
to the typical lower bound ST/T > 0.005. We have typi- 
cally used ST/T « 0.01, which is not far from the lower 
bound. For some selected temperatures, we have explic- 
itly checked that linear response is satisfied by comparing 
results for 2ST, ST and ST/2. 

It might be worth recalling that the value of XT{t) 
does not depend whether one works in NVE or NVT, 
since ensemble equivalence obviously holds for this lo- 
cal observable [91|. Much less trivial is the numerical 
finding that xt(0 is also found to be the same for New- 
tonian, Brownian and Monte Carlo dynamics for times 
pertaining to the structural relaxation. This directly fol- 
lows from the non-trivial numerical observation that the 
average structural relaxation dynamics of the binary LJ 
system has no dependence upon its microscopic dynam- 
ics, apart from an overall time rescaling. On the other 
hand, the short-time dynamics is different in the three 
cases. Our findings then confirm for Brownian dynam- 
ics, and extend for Monte Carlo dynamics [4l[, the results 
of Refs. [13, [Zil about the independence of the average 
glassy dynamics upon the microscopic dynamics. We will 
see below that clear differences emerge at the level of the 
dynamic fluctuations. 
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FIG. 4: Peak amplitude of various dynamic susceptibilities in 
the binary LJ mixture obtained from the A particles dynamics 
(top) and the BKS model for silica from the Si ions dynamics 
(bottom). Open triangles in the LJ system represent Xa^^ 
measured in a smaller system with A'' = 256 instead of the 
= 1000 used everywhere else in the paper. In both cases, 
T^Xt/cv is smaller than Xi^^ high temperature, but in- 



creases faster and becomes eventually the dominant contribu- 
tion to Xa^'^ in the relevant low temperature glassy regime. 
Note that the crossing occurs much earlier for BKS. 



C. Amplitude of the dynamic fluctuations 



In this paper, wc restrict our analysis of the dynamic 
susceptibilities to the amplitude of the peaks observed in 
Fig. [TJ meaning that we study dynamic fluctuations on 
a timescale i « Tq . The time dependence of the fluctua- 
tions are studied in the companion paper, Ref. (ssj . 



Furthermore, as discussed in Sec. Ill E 41 the contribu- 
tion to x*i due to density fluctuations in Eq. (PT|) is signif- 
icantly less than the one corresponding to energy fluctu- 
ations for most molecular liquids. Therefore, we will ne- 
glect the role of density fluctuations in the following and 
focus only on xl'^"^ since we expect that Xa^^ — Xa^^ ■ 
As a more quantitative check, we have used the data in 
Ref. [1^ to estimate that the contribution of density fluc- 
tuations to X4 is about 10 times smaller than the tem- 
perature contribution for the LJ system at po = 1.2. 



1. Ensemble dependence of dynamical correlations 

Our results are summarized in Fig. 31 where we present 
our numerical data for T\xt\/ y/cy, Xa^'^ i T'^Xt/^v 
the sum = xf^^ + T'^Xt/cv, all quantities ob- 

tained from Newtonian dynamics simulations of both the 
LJ and BKS models. Recall that we define cy in units of 
fcs throughout the paper. When temperature decreases, 
all peaks shift to larger times and track the a-relaxation. 
Simultaneously, their height increases, revealing increas- 
ingly stronger dynamic correlations as the glass transi- 
tion is approached. 

The main observation from the data displayed in Fig. 31 
already announced in Ref (37| . is that in both LJ and 
BKS systems the term T^Xt/cv while being small, 
~ O(10~^), above the onset temperature of slow dynam- 
ics, grows much faster than Xi^^ when the glassy regime 
is entered. As a consequence, there exists a temperature 
below which the temperature derivative contribution to 
the four-point susceptibility Xi^'^ dominates over that 
of Xa^^ ^ or is at least comparable. This crossover is lo- 
cated at T w 0.45 in the LJ system, T « 4500 K for BKS 
silica. The conclusion that T'^x^jcy becomes larger than 
y^^^ at low temperatures holds for both strong and 
fragile glass-formers, but for different reasons. In the LJ 
systems xt increases very fast because timescales grow 
in a super- Arrhcnius manner, which makes the temper- 
ature derivative larger and larger, while y^^^ saturates 
at low T . In the BKS system, although the tempera- 
ture derivative is not very large because of the simple 
Arrheni us g rowth of relaxation timescales, y^^^ is even 
smaller |77| . i.e. much smaller than in the fragile LJ 
system. It is interesting to note that the common value 
of Xa^^ and T^Xr/cy when they cross is substantially 
larger for the LJ system 10) than for BKS ('~-' 1). It 
would be interesting to see, more generally, how y^^^ 
and fragility are correlated. 

It is important to remark that finite size effects could 
play a role in the present study: when measured in a 
system which is too small, dynamic fluctuations arc un- 
derestimated [zll. Therefore it could be that using too 
small a system we have underestimated yJl'^^ , and there- 
fore observed a fictitious saturation of the inequality (|24p . 
To investigate this possibility, we have included in Fig. [4] 
data for y^^^ obtained in a system comprising about 
4 times less particles, TV = 256, with essentially similar 
results. We have checked that also the average dynamics 
is unchanged when N = 256, so that xt is not affected 
by finite size effects either for the range of parameters 
chosen. We are therefore confident that the main con- 
clusion drawn from Fig. HI is not an artifact due to finite 
size effects. 

We can therefore safely conclude that T^Xr/cy is an 
excellent approximation to Xa^^ f^i' relaxation times 
larger than k, 10"* in the LJ system, Tq, « 10 ps in 
BKS silica. Our results indicate that this becomes an 
even better approximation as temperature is lowered, at 
least in the numerically accessible regime. As reported in 
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FIG. 5: Amplitude of four-point susceptibilities Xi{'''a) ob- 
tained from the A particles dynamics in the LJ system for 
Newtonian canonical {X4^'^)) microcanonical {x4^^) dy- 
namics and stochastic Monte Carlo {xa'^) Brownian 
(x?^) dynamics. Stochastic dynamics measurements follow 
the results obtained from microcanonical Newtonian dynam- 
ics, while the amplitudes obtained in the canonical ensem- 
ble for Newtonian dynamics are much larger, as predicted in 
Secimi 
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FIG. 6: Four-point susceptibilities at X4(^) at T = 0.45 ob- 
tained from the A particles dynamics in the LJ system for 
Newtonian canonical (shown as a thicker line), microcanoni- 
cal dynamics and stochastic Monte Carlo and Brownian dy- 
namics as a function of a rescaled time chosen so that all X4's 
overlap near the alpha relation. We chose t = t for NVE New- 
tonian dynamics, t = t/24 for Brownian dynamics, t = t/100 
for Monte Carlo dynamics. The Newtonian Xa^'^ i't) is larger 
than the others, which are all nearly identical in both beta 
and alpha regimes. 



Ref. [131, this suggests a direct experimental determina- 
tion of X4 close to the glass transition temperature, Tg. 
Our data indicate however, that care must be taken when 
analyzing the first few decades of the dynamical slowing 
down where all terms contribute differently t o yj^ ^^, and 
have different temperature dependences [s^, [H^ [l^l . We 
now show that despite their different temperature behav- 
ior, xl'^"^, Xi^^ ^nd XT contain the same physics, as 
predicted theoretically in previous sections. 



2. Dynamics dependence of dynamical correlations 

We conclude this section with a discussion of the data 
for dynamic fluctuations obtained through our stochas- 
tic simulations. The temperature evolution of the dy- 
namic susceptibilities X'^i'^a) obtained with Monte Carlo 
and Brownian dynamics are shown in Fig. [5] where it is 
compared to the data obtained in both canonical and mi- 
crocanonical ensembles with Newtonian dynamics. Our 
data unambiguously show that dynamic fluctuations with 
stochastic dynamics are different from the ones obtained 
with Newtonian dynamics in the NVT ensemble. They 
are however very similar to the microcanonical ones. This 
result is not immediately intuitive because one could have 
imagined that stochastic simulations arc a good approx- 
imation to the dynamics of liquids in the canonical en- 
semble. However, we have shown in Sec. IIIII that this 
naive expectation is in fact incorrect. The absence of 
the energy conservation in the stochastic dynamics (MD 
or BD) removes the contribution of the "squared par- 
quets" , which corresponds to the enhancement of dy- 
namic fluctuations due to energy fluctuations, and leads 



to Xi ^ Xi ^ xi ■ This is in excellent agreement 
with our numerical data. 

Another confirmation of our theoretical expectations 
is presented in Fig. [6] in which we show the time de- 
pendence of X4 for NVE Newtonian dynamics, NVT 
Brownian, Monte-Carlo and Newtonian dynamics. The 
first three curves are essentially identical apart at micro- 
scopic times, whereas the last one is clearly larger. This 
dependence on the microscopic dynamics is a general re- 
sult obtained from the previous diagrammatic discussion. 

A further crucial prediction of our diagrammatic anal- 
ysis is that XT and X4^^ should have a similar critical 
scaling in temperature and time. This is again a general 
result if the three-leg vertex docs not introduce any addi- 
tional singular behavior. In fact, as discussed in the pre- 
vious sections, xt consists of a parquet diagram closed 
by a three-leg vertex whereas X4^^ is given by single 
parquet diagrams. In Fig. 3] we confirm numerically that 
the peaks of xt and Xi^^ scale in the same way with 
the temperature, both in the LJ and BKS systems. This 
similarity should in fact extend to the whole time depen- 
dence but the results arc somewhat less satisfactory, as 
discussed in the companion paper [ssj . 

For BKS we do not have numerical results for Brow- 
nian Dynamics for reasons mentioned above. However, 
our preliminary results from Monte Carlo simulations of a 
slightly modified version of the BKS potential [zl] agree 
with the conclusions drawn from the LJ data, that is, 
X*^*^ seems to follow more closely Xa'^^ ^ as in Fig. [51 
with similar time dependences for the dynamic suscepti- 
bilities, as in Fig. [HI 
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D. Spatial correlations 

We now discuss the spatial correlations associated 
with the global fluctuations measured through xrit) and 
X4{t)- To this end, wc define the local fluctuations of the 
dynamics through the spatial fluctuations of the instan- 
taneous value of the self intermediate scattering function, 

6fk{^,t) = ^,5(x-r,(0)) [cos[k. (r,(t) - r,(0))] - Fs{k,t)] . 0. 

(65) 

In the following, we will drop the k dependence of the 
dynamic structure factors to simplify notations. Local 
fluctuations of the energy at time t are defined as usual, 



5e(x,t) =^<5(x-r,(t))[e,(t) 



(66) 



where ei{t) — mvf{t)/2 + ^^'^ instan- 

taneous value of the energy of particle i, and e = 
{N~^ J2i s-i) is the average energy per particle. 

Spontaneous fluctuations of the dynamics can be de- 
tected through the "four-point" dynamic structure fac- 
tor, 

54(q,t) = ^(<5/(q,<)5/(-q,i)), (67) 

while correlation between dynamics and energy are quan- 
tified by the three-point function, 



5T(q,t) = ^(<5/(q,<)<5e(-q,t 



o)>. 



(68) 



In Eqs. ([6711681) . '^/(q^*) and 5e{(\,t) denote the Fourier 
transforms with respect to x of Sf{jc,t) and (5e(x, t), re- 
spectively. We will show data for fixed |k|, as for the 
dynamic susceptibilities above. In our numerical simu- 
lations we have also performed a circular averaging over 
wavevectors of fixed moduli |q|, although the relative ori- 
entations of q and k plays a role [l^, ■ 

It should be remarked that the spatial correlations 
quantified through Eqs. (|67l I68p can be measured in any 
statistical ensemble, because they are local quantities not 
sensitive to far away boundary conditions. Therefore, 
their g — > limits is related to the dynamic susceptibili- 
ties measured in the ensemble where all conserved quan- 
tities fluctuate. 

We present our numerical results for the temperature 
dependence of four-point and three-point structure fac- 
tors in Fig. [71 Similar four-point dynamic structure fac- 
tor have been discussed before [3 [l^, [H, [13, [sl], [H, 
m, [13 . They present at low q a peak whose height in- 
creases while the peak position shifts to lower q when T 
decreases. This peak is unrelated to static density fluctu- 
ations which are small and featureless in this regime [t^ . 
This growing peak is direct evidence of a growing dy- 
namic Icngthscalc, ^4(T), associated to dynamic het- 
erogeneity as temperature is decreased. The dynamic 
lengthscalc ^4 should then be extracted from these data 
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FIG. 7: Top: Four-point dynamic structure factors from 
Eqs. ^7} for Monte Carlo (MC) and Newtonian (N) dynam- 
ics, and tfiree-point structure factor, Eq. (|68p . for Newtonian 
dynamics. For comparison, we sfiow the power law as 
a dashed line. Note that St is a negative quantity, so we 
present its absolute value. St and Si''^ have been vertically 
shifted for graphical convenience. Bottom: Rescaled dynamic 
structure factor for Newtonian dynamics using Eq. (|69p with 
a = 2.4 for ^4 and a = 3.5 for St- The same dynamic length- 
scale ^4 — S^T = (. is used in both cases, and the temperature 
evolution of ^ is shown in the inset. 



by fitting the g-dependence of S4{q,t) to a specific form. 
An Ornstein-Zernike form has often been used [H, HI] , 
and we have presented its 1/g^ large q behavior in Fig. [71 
Since our primary aim is to measure dynamic suscepti- 
bilities on a wide range of temperatures, wc have used 
a relatively small number of particles, N = 1000. At 
density po = 1.2, the largest distance we can access in 
spatial correlators is L/2 « 5, which makes an absolute 
determination of ^4 somewhat ambiguous. Similarly the 
range of wavevectors shown in Fig. [71 is too small to as- 
sign a precise value even to the exponent characterizing 
the large q behavior of 5*4(9, t) ~ l/?"- Our data is com- 
patible with a value a « 2.4. To extract ^4 we therefore 
fix a = 2.4 and determine ^4 by assuming the following 
scaling behavior [33 |. 



(69) 



usm 



ig 5*4(0, i) and ^4 as free parameters. The results of 
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such an analysis are shown in the bottom panel of Fig. [T] 
This procedure leads to values for ^4 which are in good 
agreement with previous determinations using different 
procedures [11]. In particular we find that a power law 

\l z 

relationship ^4 ~ Tq with z « 4.5 describes our data 
well, as reported in Ref. [l^ for this system. 

Since dynamic structure factors probe local spatial cor- 
relations they do not depend on the statistical ensemble 
chosen for their calculation, at least in the thermody- 
namic limit. However, as predicted in Sec. IIIII dynamic 
correlations are expected to retain a dependence on the 
microscopic dynamics of the particles, our prediction be- 
ing that correlations should be stronger for Newtonian 
dynamics than for stochastic dynamics. This prediction 
is directly confirmed in Fig. [7| where we show 84:'^'-' {q,t) 
obtained from our Monte Carlo simulations. Clearly the 
temperature evolution of 5*4^'^ is slower than that of 
S^^, in agreement with the slower temperature evolu- 
tion of already observed in Fig. [S] 

An important new result contained in Fig. [7] is the 
presence and development of a similar low-q peak in the 
three-point structure factor ST{q,t). Note that, as for 
Xrit), we find that 5't((7, i) is a negative quantity. This 
means that a local positive fluctuation of the energy is 
correlated to a local negative fluctuation of the two-time 
dynamics, i.e. to a locally faster than average dynam- 
ics. Therefore the (negative) peak in ST{q,t) is a direct 
microscopic demonstration that dynamic heterogeneity 
is strongly correlated to the fluctuations of at least one 
local structural quantity, namely the energy [sj. When 
temperature decreases, the height of the peak in |S't((J, t)\ 
increases and it shrinks towards lower q. This is again the 
sign of the presence of a second growing dynamic length- 
scale, (,T, which reflects the extent of the spatial correla- 
tions between energy and dynamical fluctuations. Again 
an absolute determination of is very hard due to sys- 
tem size limitations. Since we expect and ^4 to carry 
equivalent physical content, we have checked that our 
data are compatible with both lengthscales being equal. 
In Fig. [3 we rescale the three-point dynamic structure 
factor using Eq. (j69p with a = 3.5, and constraining 
?r = ^4- The scaling is of similar quality, sec the bot- 
tom panel in Fig.[7l Clearly, the nontrivial g-dependence 
of ST{q,t) with a scaling collapse of reasonable quahty 
and a length scale consistent with that extracted from 
S4{q,t) is a strong indicator that the integrated suscep- 
tibility Xt grows as a result of a unique growing length 
scale characteristic of dynamic heterogeneity. A further 
numerical confirmation of the fact that the growth of 
the susceptibilities Xii^) and xt(0 cannot be attributed 
to an increase in the strength of the correlations rather 
than their range stems from the direct measurement of 
(74 (r, t) and grir, t), the Fourier transforms of 6*4(9, t) and 
ST{q,t). The large distance decay of both functions can 
be well fitted, within the statistical noise, by an exponen- 
tial form with a growing dynamic lengthscale but a 
temperature independent strength. Attributing all of the 
temperature dependence of the susceptibilities to a grow- 



ing amplitude leads to poor fits of the spatial correlators. 
This indicates that a scenario whereby the growth of xt 
can be ascribed to the growth of a prefactor with no grow- 
ing length scale characteristic of dynamic heterogeneity 
can be ruled out, at least for the LJ case. Collectively, 
these findings indicate that the bound for Xi, as first dis- 
cussed in [37[, correctly estimates a correlation volume 
associated with dynamic heterogeneity. 

We have carried out a similar analysis for the BKS 
model of silica [t^. Here, the analysis is far more difficult 
for several reasons. First, the system is harder to simu- 
late than the LJ system due to the long-ranged nature of 
the interactions. Second, strong features associated with 
static structure make a resolution of the low-g behavior 
in STiq,t) somewhat more challenging in this system. 
Lastly, the overall scale of dynamical fluctuations at the 
lowest temperatures studied are much smaller than in 
the LJ system (see Fig. 2]). Regardless, we do find re- 
sults consistent with a scaling scenario for ST{q,t)^ and, 
as we will see in the following paper, the growth of xt 
tracks that of Xi^'^- These facts give support to the no- 
tion that the scenario for the BKS model of silica is the 
same as for the LJ system although the direct supporting 
evidence for this is, at this stage, not quite as strong. 

The local correlation between energy fluctuations dy- 
namic heterogeneity is broadly consistent with several 
theoretical predictions; see the companion paper [sst for 
further discussion. As mentioned in Sec. IIIII the equality 
between and ^4 is a natural prediction, in particular 
close to the MCT transition. This is also very natural 
from the point of view of kinetically constrained mod- 
els j44j . Spin facilitated models, in particular, postulate 
such a correlation through the concept of dynamic fa- 
cilitation: mobile sites carry positive energy fluctuations 
and through activated diffusion trigger the relaxation of 
neighboring sites [sij . In this picture a localized energy 
fluctuation affects the dynamics of a large nearby region 
so that there is not a one-to-one correspondence between 
slow and low-energy sites. There is therefore no contra- 
diction between our results and the lack of correlation 
between "dynamic propensity" and local potential en- 
ergy recently reported in Ref. (s^ . They qualitatively 
agree, however, with recent numerical results obtained 
for water where a correlation between "dynamic" and 
"energetic" propensities is reported [s^l- A recent work 
[ssj has also suggested a relation between energy fluctu- 
ations and finite size effects leading to a growing length 
at low temperature. 



E. Summary 

In this section we have discussed in detail the results of 
molecular dynamics simulations of a strong and a frag- 
ile glass-forming liquid. Our main contribution is the 
simultaneous measurement of spontaneous and induced 
dynamic fluctuations, and the quantitative conflrmation 
in two realistic liquids of the central claim announced 
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in Ref. [37[: it is possible to obtain a quantitative esti- 
mate of the amplitude of dynamic fluctuations in super- 
cooled liquids through the measurement of the quantity 
T^Xt I yfcy ■, which (once squared) gives the major contri- 
bution to Xa^'^ : and hence to X4^^^ i^i the low temper- 
ature regime, and is proportional to X4^^- see Fig. S) 

We have directly measured in the Lennard-Jones sys- 
tem three- and four-point dynamic structure factors that 
display slightly different wavevcctor dependences but 
lead nevertheless to consistent quantitative estimates of 
a dynamic correlation lengthscalc. compatible with that 
obtained from and xt- This last result is very impor- 
tant since this is direct confirmation that an experimen- 
tal estimate of a dynamic lengthscale, as performed in 
Ref. [13, is meaningful. Finally, we have found that, as 
predicted theoretically, global four-point dynamic corre- 
lations corresponding to spontaneous fluctuations of two- 
time correlators are strongly dependent on the micro- 
scopic dynamics, at variance with usual two-point corre- 
lations. 



VI. PERSPECTIVES AND CONCLUSIONS 

We conclude this rather long article, to be followed 
by a companion paper [ssj . with brief comments only. 
Four-point correlators were originally introduced to de- 
fine the lengthscale of dynamical heterogeneities in glass 
formers. Our results, in that respect, are double-sided. 
We showed that global four-point functions, correspond- 
ing to the fluctuations of intensive dynamical correlators, 
not only depend on the statistical ensemble but also on 
the choice of dynamics. The dependence on the statis- 
tical ensemble is useful to obtain lower bounds for ex- 
perimentally relevant situations. However, on a more 
general ground, these dependences unveil that four-point 
correlators are more complicated than what was orig- 
inally thought and their quantitative interpretation is 
somewhat flummoxed. Instead, we found that dynamical 
response functions, proxied by the temperature or den- 
sity derivatives of two-time correlators provide a more 
clear and direct probe of genuine collective dynamical ef- 
fects. 

We have given strong theoretical and numerical evidence 
for the most important claim made in this paper: the ex- 
istence of unique dynamical lengthscale ^ governing the 
growth of all the relevant dynamical susceptibilities, in- 
dependently of dynamics (and of course ensemble!). This 



result can be proved within the MOT of glasses, as we 
elaborate further in the companion paper (38[ , but is ex- 
pected more generally as soon as ^ becomes somewhat 
large compared to the inter-atomic spacing. Our numer- 
ical results show that this is true both in the fragile LJ 
system and in the strong BKS system: all dynamical sus- 
ceptibilities {xt, Xi^^) behave similarly, at least in the 
weakly supercooled region accessible to numerical simu- 
lations. 

One rather striking result of our analysis is that even 
Arrhcnius dynamics in Newtonian systems must involve 
some amount of dynamical correlations. This is con- 
firmed by our numerical simulations on the BKS system, 
but the result holds more generally. Even a dilute assem- 
bly of Arrheniusly relaxing entities, e.g. two-level sys- 
tems, should develop non-trivial dynamical correlations 
at sufficiently low temperatures, provided they interact 
with the same Newtonian thermal bath. This is obvi- 
ously the case for a strong glass-formers where a particle 
is both a relaxing entity and part of its neighbors' bath. 

As for the perspective for the future, we hope that 
our work will trigger more experimental and numerical 
investigations of supercooled liquids and jamming sys- 
tems, extending our results both from a quantitative and 
a qualitative point of view (45j . In particular, the dis- 
tinction between dynamical correlations (explored here) 
and cooperativity, if any, should be clarified. The rela- 
tion between the two notions might be very different in 
strong and fragile systems, and the distinction between 
the MCT and the deeply supercooled regimes might also 
be relevant. Is ^ as defined in the present paper related to 
the Adam-Gibbs or the mosaic length scale [13, [sBl? In 
this respect, the full understanding of deceivingly simple 
Arrhcnius systems should be of great help. 
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